Impact of a global quadratic potential on galactic rotation curves 

Philip D. Mannheim and James G. O'Brien 
Department of Physics 
University of Connecticut 
Starrs, CT 06269, USA 
^1 philip.mannheim@uconn.edu, 

>• ■ obrienOphys.uconn.edu 

^ ■ (Dated: November 21, 2010) 

CN : Abstract 

^^ ■ We have made a conformal gravity fit to an available sample of 110 spiral galaxies, and report 

^^ ! here on the 20 of those galaxies whose rotation curve data points extend the furthest from galactic 






o 



X 



centers. We identify the impact on the 20 galaxy data set of a universal de Sitter-like potential term 



O ' V{r) = — Kc^r^/2 that is induced by inhomogeneities in the cosmic background. This quadratic 

t^ . term accompanies a universal linear potential term V{r) = 7oc^r/2 that is associated with the 

cosmic background itself. We find that when these two potential terms are taken in conjunction with 

<N 

^ ' the contribution generated by the local luminous matter within the galaxies, the conformal theory is 

O' 

C^^ ■ able to account for the rotation curve systematics that is observed in the entire 110 galaxy sample, 

Q\ 

^^ \ without the need for any dark matter whatsoever. With the two universal coefficients being found to 

|> 

^^ ' be of global magnitude, viz. k = 9.54 x 10 ^'^ cm ^ and 70 = 3.06 x 10 ^"^cm ^, our study suggests 



that invoking the presence of dark matter may be nothing more than an attempt to describe global 
effects in purely local galactic terms. With the quadratic potential term having negative sign. 



H \ galaxies are only able to support bound orbits up to distances of order 7o/k = 3.21 x 10^'^ cm, with 

global physics thus imposing a natural limit on the size of galaxies. 



I. INTRODUCTION 

At the present time it is widely believed that on scales much larger than solar-system- 
sized ones astrophysical and cosmological phenomena are controlled by dark matter and 
dark energy, with luminous matter being only a minor contributor. However, given the lack 
to date of either direct detection of dark matter particles or of a solution to the cosmological 
constant problem, a few authors (see e.g. [1] for a recent review) have ventured to suggest 
that the standard dark matter/dark energy picture may be incorrect, and that one instead 
needs to modify the standard Newton-Einstein gravitational theory that leads to that picture 
in the first place. In this paper we study one specific alternative to Einstein gravity that 
has been advanced, namely conformal gravity. We report here on the results of a conformal 
gravity study of the instructive 20 largest of a full sample of 110 galaxies, all of whose 
rotation curves we have been able to fit without the need for any dark matter at all. 

In seeking an alternative to Einstein gravity that is to address both the dark mat- 
ter and dark energy problems, our strategy is to seek some alternate, equally metric- 
based theory of gravity that possesses all of the general coordinate invariance and equiv- 
alence principle structure of Einstein gravity, that yields a geometry that is described 
by the Ricci-fiat Schwarzschild metric on solar-system-sized distance scales while depart- 
ing from it on larger ones, and that has a symmetry that can control the cosmological 
constant A. All of these criteria are met in the conformal gravity theory (see e.g. [l|) 
that was first developed by Weyl. Specifically, as well as coordinate invariance, in ad- 
dition one requires that the action be left invariant under local conformal transforma- 
tions of the form g^,y{x) — )■ e'^"^^^g^^{x) with arbitrary local phase a{x). Given this re- 
quirement, the gravitational action is then uniquely prescribed to be of the form Jw = 

-«,/rf4x(-(?)l/2CA^,,C^'^- = -2«,/rf4x(-(7)l/2[i?^,i?'^--(l/3)(i?°j2] ^l^g^g c^f..n -^ 

the conformal Weyl tensor and ag is a dimensionless gravitational coupling constant. 

With the conformal symmetry forbidding the presence of any fundamental A term in Jw, 
conformal gravity has a control on A that is not possessed by Einstein gravity; and through 
this control conformal gravity is then able to solve the cosmological constant problem |2|, yj. 
In addition, the conformal gravity equations of motion are given by [l| 



AagW^"' = Aag 



r\/^^l\uK /^flXuK, p 



X:k-'-' -f^AK 



4ag 



W"(2)-^W^WI =^^^ 



Wl"^. = 2g^^{R\y^, - 2{R\y^-^^ - 2R\R^^ + -g^'^iR^f 



(1) — -^y V-'i' a; ■l3~^\-'^a) - ^-n. Q,-". -r -I 



fl 



with Schwarzschild thus being a vacuum solution to conformal gravity, just as required J4|. 



II. UNIVERSAL POTENTIALS FROM THE REST OF THE UNIVERSE 

Given its structure, W^'^ could potentially vanish even if the geometry is not Ricci flat, 

and the conformal theory could thus have non-Schwarzschild solutions as well. To identify 

any such solutions, Mannheim and Kazanas solved for the metric in a vacuum region exterior 

to a static, spherically symmetric source of radius tq, to find 51] that the exact, all-order line 

element is given by ds^ = —B{r)dt^ + dr"^ / B{r) + r'^dQ2, with 

2/3 

B(r > ro) = 1 + 7r - kr'^. (2) 

r 

In the jr and —kr"^ terms we see that the conformal gravity metric departs from the B{r > 
ro) = 1 — 2/3/r Schwarzschild metric only at large r, just as we want. 

I. seeking to relate the various constants m O to properfes of the e,rergy-,nonrentum 
tensor T^,^ of the source, Mannheim and Kazanas found [6] that in the static, spherically 
symmetric case the quantity (3/-B(r)) {W^q — W^^) evaluates exactly to V^B{r)] and that, 
in terms of the general source function /(r) = (3/4a;gi?(r))(T''Q— T^^), the exact fourth-order 
equation of motion of the conformal theory reduced to the remarkably simple 



V^B 



d^ 4 d^ 
dr^ r dr^ 



B{t) = f{r), (3) 



without any approximation whatsoever. Since V'^(r^) vanishes identically everywhere while 
V^(l/r) and V^(r) evaluate to delta functions and their derivatives, we see that of the 
constants given in (|2]), only /3 and 7, but not k, can be associated with properties of a local 
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source of radius ro; with the matching of the interior and exterior metrics then yielding |6| 

1 pro 1 pro 

7 = --/ rfrVV(r'), 2/3 = -/ drV'fir'). (4) 
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Thus despite the presence of the — fcr^ term in the exterior vacuum solution in ([2]), the 
above analysis provides no specific basis for considering it further, as it is associated with 
the trivial solution to V^-B = 0, to thereby be devoid of dynamical content. 



In conformal gravity a local gravitational source generates a gravitational potential 

per unit solar mass, with /3* being given by the familiar MqG/c^ = 1.48 x 10^ cm, and 
with the numerical value of the solar 7* needing to be determined by data fitting. In the 
theory the visible local material in a given galaxy would generate a net local gravitational 
potential VlocI'^) given by integrating V*{r) over the visible galactic mass distribution. In 
disk galaxies luminous matter is typically distributed with a surface brightness S(i?) = 
Sqc"'^/'^'' with scale length Rq and total luminosity L = 27iT.oRq, with most of the surface 
brightness being contained in the R < 4Rq or so optical disk region. For a galactic mass to 
light ratio M/L, one can define the total number of solar mass units A^* in the galaxy via 
{M/L)L = M = N*Mq. Then, on integrating V*{r) over this visible matter distribution, 
one obtains [1| the net local luminous contribution 
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for the centripetal accelerations of particles orbiting in the plane of the galactic disk. 

However, unlike the situation that obtains in standard second-order gravity, one cannot 
simply use ([6]) as is to fit galactic rotation curve data, as there are two additional global 
effects coming from the rest of the material in the universe that need to be taken into 
consideration as well, one associated with the homogeneous cosmological background and 
the other with the inhomogeneities in it. As regards first the effect of inhomogeneities, we 
recall for the standard second order Poisson equation V^(j){r) = g{r), the force associated 
with a general static, spherically symmetric source g{r) is given by 



2 dr'r"g{r'). (7) 

ar r^ Jo 

As such, the import of (JTj) is that even though g{r) could continue globally all the way 
to infinity, the force at any radial point r is determined only by the material in the local 
< r' < r region. In this sense Newtonian gravity is local, since to explain a gravitational 
effect in some local region one only needs to consider the material in that region. Thus in 



Newtonian gravity, if one wishes to explain the behavior of galactic rotation curves through 
the use of dark matter, one must locate the dark matter where the problem is and not 
elsewhere. Since the discrepancy problem in galaxies occurs primarily in the region beyond 
the optical disk, one must thus locate galactic dark matter in precisely the region in galaxies 
where there is little or no visible matter. 

However, this local character to Newtonian gravity is not a generic property of any 
gravitational potential. In particular for the fourth-order Poisson equation V^0(r) = h{r) = 
/(r)c^/2 of interest to conformal gravity, the potential and the force evaluate to 

<^(r) = — / dr'r''^hir') / dr'r'%ir') - - / dr'r'^hir') / dr'r'hir'), 
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'^'^^'"^ - -- r dr'r'^h{r') + ^ r dr'r'^h{r')-- I dr'r'h{r'), (8) 
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so that this time we do find a global contribution to the force coming from material that is 
beyond the radial point of interest. Hence in conformal gravity one cannot ignore the rest 
of the universe, with a test particle in orbit in a galaxy being able to sample both the local 
field due to the matter in the galaxy and the global field due to the rest of the universe. 

In the presence of inhomogeneities W'^''' does not vanish, as the very presence of a localized 
source prevents a geometry from being conformal to flat, with inhomogeneities in the universe 
thus leading to integrals in ([H]) that can extend to very large distances. However, this is not 
the only global effect that we need to take into consideration, as one can also add on to 
(j8]) any terms that would cause W'^'^ to vanish, provided they make it do so non-trivially. 
Since the cosmological Robertson- Walker (RW) metric is homogeneous and isotropic, it is 
conformal to fiat, and thus its geometry obeys W^" = 0. For the cosmological background 
the vanishing of W^" entails that conformal cosmology be described by T^"^ = 0. As discussed 
in [l| the equation T^'^ = can be satisfied non-trivially, and leads to a topologically open 
RW cosmology, with its contribution to W'^'^ then vanishing non-trivially, just as desired. 

Since cosmology is written in comoving Hubble flow coordinates while rotation curves are 
measured in galactic rest frames, to ascertain the impact of cosmology on rotation curves one 
needs to transform the RW metric to static coordinates. As noted in [Sj], the transformation 

4r f 

P = 2(l + 7or-fcr^)V^ + 2 + 7or ' " = J ^'"^^'^ ^'^ 

effects the metric transformation 

dv 

-(1 + 7or - kr^ydt^ + J— — - + r^dQ^ = 

(1 -I- 7or — /cr^) 
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Since an RW geometry is conformally flat and since it remains so under a conformal trans- 
formation, we see tliat wlien written in a static coordinate system, a comoving conformal 
cosmology with 3-space spatial curvature K looks just like a static metric with universal 
linear and quadratic terms with coefficients that obey K = — 7o/4 — k. However, since there 
was only one spatial scale in the RW metric (viz. K), its decomposition into two static 
coordinate system scales (70 and k) was artificial, and so in [7] the k term was dropped. 
Then, without the k term, we see that in the rest frame of a comoving galaxy (i.e. one with 
no peculiar velocity with respect to the Hubble fiow), a topologically open RW cosmology 
would look just like a universal linear potential with cosmological strength 7o/2 = (—KY^'^. 
In the conformal theory then we recognize not one but two linear potential terms, a local 
A^*7* dependent one associated with the matter within a galaxy and a global cosmological 
one 7oc^r/2 associated with cosmological background. Thus in the weak gravity limit one 
can add the two potentials and replace ([6]) by [3] 
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In [7J (ITT]) was used to fit the galactic rotation curve data of a sample of 11 galaxies (of 
which only NGC 2841 and NGC 3198 are in the sample considered here), and good fits were 
found, with the two universal linear potential parameters being fixed to the values 

7* = 5.42 X 10-^^cm-\ 70 = 3.06 x lO'^^cm'^ (12) 

The value obtained for 7* entails that the linear potential of the Sun is so small that there 
are no modifications to standard solar system phenomenology, with the values obtained for 
A^*7* and 70 being such that one has to go to galactic scales before their effects can become 
as big as the Newtonian contribution. The value obtained for 70 shows that it is indeed of 
cosmological magnitude. In the fitting to the 110 galaxy sample flT2|) does not change. 

However, as we had noted above, there is a contribution due to inhomogeneities in the 
cosmic background that we need to include too. These inhomogeneities would typically be 
clusters and superclusters and would be associated with distance scales between 1 Mpc and 
100 Mpc or so. Without knowing anything other than that about them, we see from (jHj) that 
for calculating potentials at galactic distance scales (viz. scales much less than cluster scales) 



the inhomogeneities would contribute constant and quadratic terms multiplied by integrals 
that are evaluated between end points that do not depend on the galaxy of interest, to thus 
be constants. Thus, again up to peculiar velocity effects, we augment ( |TT]) to 

% = % + ^ - ^c'R, (13) 

with asymptotic limit 

%-^ + ^.^-.c^«. (14) 

It is thus (TT3|) with its universal k that we must use for fitting galactic rotation curves, and 
in making such fits the only parameter that can vary from one galaxy to the next is the 
galactic disk mass to light ratio as embodied in A^*. Our fits are thus highly constrained, 
one parameter per galaxy, fits (the fits also include the effect of HI gas, but for the gas the 
mass is known), with everything else being universal, and no dark matter being assumed. 

III. DATA FITTING 

We recall that in [7i| successful rotation curve fitting to an 11 galaxy sample was obtained 
using f lTTj) . and one would thus initially anticipate that even if the —kc^R term were to 
be present in principle, in practice it would be too small to have any effect. However, the 
sample we have studied now is altogether larger (110 galaxies) and it contains some very 

instructive galaxies whose data points extend to larger distances from galactic centers than 

p 

had been the case for the 11 galaxy sample studied in [7]. It is through fitting these highly 
extended galaxies that we are able to uncover a role for the —nc^R term and extract a value 
for K given by k = 9.54 x 10~^^ cm~^. And in the fitting to the full 110 galaxy sample 



B 



to be reported elsewhere [8| (a varied sample of galaxies that includes high (HSB) and low 
(LSB) surface brightness galaxies and dwarfs) we are able to confirm that even with this 
now fixed value for n, ( IT3|) fully accounts for the data . With k being found to be of order 
1/(100 Mpc)^, it is indeed an inhomogeneous rather than a Hubble distance scale. 

In Fig. (1) we present our fits to the 20 galaxy sample with the relevant parameters being 
listed in Table (1). In Fig. (1) the rotational velocities and errors (in km sec~^) are plotted 
as a function of radial distance (in kpc) . For each galaxy we exhibit the contribution due 
to the luminous Newtonian term alone (dashed curve), the contribution from the two linear 



terms alone (dot dashed curve), the contribution from the two hnear terms and the quadratic 
terms combined (dotted curve), with the full curve showing the total contribution. As we 
see, without any need for dark matter, our fitting captures the essence of the data. Because 
the data go out to much further distances than had been the case for the sample studied 
in |7|, the data are now sensitive to the rise in velocity associated with the linear potential 
terms, and it is here that the quadratic term acts to actually arrest the rise altogether 
(dotted curve) and cause all rotation velocities to ultimately fall. Moreover, since f ^ cannot 
be negative, beyond a distance R of order 7o/k = 3.21 x 10^^ cm or so there could no longer 
be any bound galactic orbits, with galaxies thus having a natural way of terminating, and 
with global physics thus imposing a natural limit on the size of galaxies. To illustrate this 
we plot the rotation velocity curve for UGC 128 over an extended range. The fits presented 
here and in p] are noteworthy since the universal 70 and n terms have no dependence on 
individual galactic properties whatsoever and yet have to work in every single case. 



It is important to appreciate that the fits provided by con 



brmal gravity (and likewise 



by other alternate theories such as the MOND |9| and MSTG lOj] theories are predictions. 
Specifically, for all these theories the only input one needs is the optical data, and the only 
free parameter is the M/ L ratio for each given galaxy, with rotation velocities then being 
determined pJJ]. As Table (1) shows, by and large the M/ L ratios that are found are all of 
order the solar Mq/Lq ratio, just as one would want. It is important to emphasize that the 
fits are predictions since dark matter fitting to galactic data works very differently. There one 
first needs to know the velocities so that one can then ascertain the needed amount of dark 
matter, i.e. in its current formulation dark matter is only a parametrization or postdiction 
of the velocity discrepancies that are observed and is not a prediction of them. Dark matter 
theory has yet to develop to the point where one is able to predict rotation velocities given a 
knowledge of the luminous distribution alone. Thus dark matter theories, and in particular 
those theories that produce dark matter halos in the early universe, are currently unable to 
make an a priori determination as to which halo is to go with which particular luminous 
matter distribution, and need to fine-tune halo parameters to luminous parameters galaxy 
by galaxy. No such shortcoming appears in conformal gravity, and if standard gravity is 
to be the correct description of gravity, then a universal formula akin to the one given in 
flT3|) would need to be derived by dark matter theory. However, since our study establishes 
that global physics does indeed influence local galactic motions, the invoking of dark matter 
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in galaxies could potentially be nothing more than an attempt to describe global effects in 
purely local galactic terms. We would like to thank Dr. J. R. Brownstein, Dr. W. J. G. de 
Blok, Dr. J. W. Moffat, and Dr. S. S. McGaugh for helpful communications, and especially 
for providing their galactic data bases. We are particularly indebted to Dr. McGaugh for 
having alerted us to the fact that a linear potential would lead to an overshoot in UGC 128. 
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TABLE I: Properties of the 20 Large Galaxy Sample 



Galaxy Type Distance Lb Ro Rust Mm M^isk (M/L)stars (i' /c -R)iast 

(Mpc) (IOIOL0) (kpc) (kpc) (IOIOM0) (10^0^0) {MqL-^) (10-30cm-i) 



NGC 3726 


HSB 


17.4 


3.34 


3.2 


31.5 


0.60 


3.82 


1.15 


3.19 


NGC 3769 


HSB 


15.5 


0.68 


1.5 


32.2 


0.41 


1.36 


1.99 


1.43 


NGC 4013 


HSB 


18.6 


2.09 


2.1 


33.1 


0.32 


5.58 


2.67 


3.14 


NGC 3521 


HSB 


12.2 


4.77 


3.3 


35.3 


1.03 


9.25 


1.94 


4.21 


NGC 2683 


HSB 


10.2 


1.88 


2.4 


36.0 


0.15 


6.03 


3.20 


2.28 


UGC 1230 


LSB 


54.1 


0.37 


4.7 


37.1 


0.65 


0.67 


1.82 


0.97 


NGC 3198 


HSB 


14.1 


3.24 


4.0 


38.6 


1.06 


3.64 


1.12 


2.09 


NGC 5371 


HSB 


35.3 


7.59 


4.4 


41.0 


0.89 


8.52 


1.44 


3.98 


NGC 2998 


HSB 


59.3 


5.19 


4.8 


41.1 


1.78 


7.16 


1.75 


3.43 


NGC 5055 


HSB 


9.2 


3.62 


2.9 


44.4 


0.76 


6.04 


1.87 


2.36 


NGC 5033 


HSB 


15.3 


3.06 


7.5 


45.6 


1.07 


0.27 


3.28 


3.16 


NGC 0801 


HSB 


63.0 


4.75 


9.5 


46.7 


1.39 


6.93 


2.37 


3.59 


NGC 5907 


HSB 


16.5 


5.40 


5.5 


48.0 


1.90 


2.49 


1.89 


3.44 


NGC 3992 


HSB 


25.6 


8.46 


5.7 


49.6 


1.94 


13.94 


1.65 


4.08 


NGC 2841 


HSB 


14.1 


4.74 


3.5 


51.6 


0.86 


19.55 


4.12 


5.83 


UGC 0128 


LSB 


64.6 


0.60 


6.9 


54.8 


0.73 


2.75 


4.60 


1.03 


NGC 5533 


HSB 


42.0 


3.17 


7.4 


56.0 


1.39 


2.00 


4.14 


3.31 


NGC 6674 


HSB 


42.0 


4.94 


7.1 


59.1 


2.18 


2.00 


2.52 


3.57 


UGC 6614 


LSB 


86.2 


2.11 


8.2 


62.7 


2.07 


9.70 


4.60 


2.39 


UGC 2885 


HSB 


80.4 


23.96 


13.3 


74.1 


3.98 


8.47 


0.72 


4.31 
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FIG. 1: Fitting to the rotational velocities of the 20 galaxy sample 
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